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A neutrino event generator developed by the Wroclaw Neutrino Group is described. The physical models 
included in the generator are discussed and illustrated with the results of simulations. The considered processes 
are quasi-elastic scattering and pion production modelled by combining the A resonance excitation and deep 
inelastic scattering. 



1. Introduction 

Many of the neutrino experiments now oper- 
ating or in preparation are those with artificial 
neutrino beams. The beam energy ranges from 
0.5 to about 30 GeV (MINOS, CNGS, T2K). The 
uniform theoretical description of neutrino inter- 
action within this region is problematic, because 
three different dynamics: quasielastic(QE), reso- 
nance excitation(RES), and deep inelastic scat- 
tering(DIS) must be taken into account. Espe- 
cially in the region of about 1 GeV they con- 
tribute with comparable strength and the prob- 
lem of the overlap between the DIS formalism and 
the resonances production model appears. 

Up to now most of the Monte Carlo codes 
use the Rein-Sehgal model of the resonances pro- 
duction. However, another approach has been 
proposed X which makes use of the postulated 
quark-hadron duality in neutrino reactions. One 
can assume that the contributions from higher 
resonances are averaged by the deep inelastic 
scattering structure functions and only the domi- 
nant A resonance has to be treated separately. In 
this approach the exclusive inelastic final states 
must be obtained from the DIS formalism even 
in the resonance region W 6 (M + m^, 2 GeV), 
where the most important contribution comes 
from pion production channels. 



2. Generator description 

Dynamical models included in the generator 
are described below. For the quasielastic scat- 



tering the usual Llewellyn Smith |2J model is 
used with BBBA05 form factors pj]. For elas- 
tic scattering new strange form factors are taken 
into account [3]. For antineutrino CC interaction 
there are three additional channels with hyperons 
A, S+, E° production!] AY" | = 1). The A reso- 
nance is described by the form factors as in 0, 
and for DIS we use parton distribution functions 
GRV94 [H] with modification proposed by jjj. In 
order to get the final state from the inclusive 
DIS cross section we use model of fragmentation 
based on the LUND model, implemented in the 
PYTHIA6 generator. It is assumed that the in- 
teraction occurs on a separate constituent of the 
nuclcon. The fragmentation procedure was de- 
scribed in 8 . In order to get good agreement 
with experimentally measured charged particle 
multiplicity [§] we have changed a few PYTHIA6 
parameters. The final values, fine-tuned for spe- 
cific procedure of choosing the interacting parton 
inside the nucleon, are : PARJ(32)=0.1 GeV, 
PARJ(33) = 0.5 GeV, PARJ(34) = 1 GeV, 
PARJ(36)= 0.3 GeV, MSTJ(17) = 3. The defi- 
nitions and default values of the parameters one 
can find in PYTHIA manual [Ej. 

We use a procedure which gives exclusive cross 
sections for each single pion production channel 
(SPP): 



da 



SPP 



da L 



(l-a(W)) 



da DIS 



^SPP 



(W)a(W) 



dW dW v v " dW 

where F SPP is the percentage of the given SPP 
channel within the overall DIS cross section. The 
function a(W) defines transition between 
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Figure 1. Cross section for neutrino scattering 
vN -> iTX. 
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Figure 4. Cross section for double-pion produc- 
tion: vn — > [i~pTr + ir~ . 




Figure 2. Cross section for antmeutrmo scatter- _ „ , . . . 
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two SPP models: mediated by the A resonance 
and the inelastic one, and accounts for the non- 
resonant background. 

3. Results 

Many of the results from the Wroclaw Neutrino 
Generator have already been presented jHj- The 
basic results for CC neutrino and antineutrino 
scattering off nucleon are shown in figs.^121)- The 
cross section for SPP channel induced by antineu- 
trino is shown in fig- EI In figs. ^ 121 and |3 only 
SPP events with W < 2 GeV are taken into ac- 
count. The cross section for two channels with 
two final pions, for which experimental data ex- 
ists are presented in figs.|l||5J and the sum of the 
cross sections for three pions production is given 
in fig. H3 The simulations results agree with the 
experimental data for single and double pion pro- 
duction, but some discrepancy occurs for triple- 
pion production. 
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